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We have measured the ratios of antiparticles to particles 
for charged pions, kaons and protons near mid-rapidity in 
central Au+Au collisions at ^/s NN = 130 GeV. For protons, 
we observe (f)/(p) = 0.60 ± 0.04(staf.) ± 0.06(syst.) in the 
transverse momentum range 0.15 < pr < 1.0 GeV/c. This 
leads to an estimate of the baryo-chemical potential of 
45 MeV, a factor of 5-6 smaller than in central Pb+Pb 
collisions at ^/s NN = 17.2 GeV. 

PACS numbers: 25.75.-q 

In this paper multiplicity ratios of antiparticles to 
particles for primary charged pions, kaons and protons 
are presented for collisions of gold nuclei at an energy 
of y/s NN — 130 GeV. The data were taken with the 
PHOBOS detector during the first run of the Relativistic 
Heavy-Ion Collider (RHIC) at Brookhaven National Lab- 
oratory. The experiments at RHIC aim at understanding 
the behavior of strongly interacting matter at high 
temperature and density. Quantum chromodynamics, 
the fundamental theory of strong interactions, predicts 
that under these conditions a new state of matter will be 
formed, the quark-gluon plasma yj] . 

One of the most intriguing results from heavy-ion col- 
lisions at lower energies was the observation that particle 
ratios for particles with production cross-sections varying 
by several orders of magnitude could be described in a 
statistical picture of particle production assuming chem- 
ical equilibrium This is particularly remarkable 
for the production rates of rare baryons and antibaryons 
containing multiple strange quarks § , which were found 
to be difficult to reproduce in microscopic hadronic 
transport models Hf7|. One of the key ingredients for 
the statistical particle production picture is the baryo- 
chemical potential ps- The particle ratios presented 
here, in particular (p)/(p), provide information about ps 
at RHIC energies. 



On a microscopic level, the (p)/(p) ratio reflects 
the interplay between the transport of baryon number 
carried by protons and neutrons in the colliding nuclei 
in momentum space, the production of quark-antiquark 
pairs and the annihilation of antiprotons in the final 
stages of the collision. The connection between baryon 
number transport and the energy loss of the incoming 
nucleons [|| is under intense theoretical discussion ||,[l0) . 

The data reported here were collected using the PHO- 
BOS two-arm magnetic spectrometer. Details of the 
setup can be found elsewhere ^TJ-413). One arm (SPECP) 
of the spectrometer was only partially equipped with 6 
layers of silicon sensors, providing tracking only in the 
field-free region close to the beampipe. The other arm 
(SPECN) had a total of 16 layers of sensors, providing 
tracking both outside and inside the 2 T field of the 
PHOBOS magnet. Particles within the geometrical 
acceptance region used in this analysis traverse 14 or 15 
of the layers. A two layer silicon detector (VTX) covering 
\rj\ < 1.5 and 25% of the azimuthal angle provided addi- 
tional information on the position of the primary collision 
vertex. In total, 94362 sensitive detector elements were 
read out, of which less than 2% were non- functional. 

The primary event trigger was provided by two sets of 
16 scintillator paddle counters covering pseudorapidities 
3 < \rj\ < 4.5. Additional information for event 
selection was obtained from two zero-degree calorimeters 
measuring spectator neutrons. Details of the event 
selection and centrality determination can be found in 
Jl4] , p^[ . Monte Carlo (MC) simulations of the apparatus 
were based on the HIJING event generator jl(| and the 
GEANT 3.21 simulation package, folding in the signal 
response for scintillator counters and silicon sensors. 

For this analysis the 12% events with the highest signal 
in the paddle counters were selected, corresponding to 
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the collisions with the highest number of participating 
nucleons N part . The average number of participants for 
the selected events was estimated as (N part ) = 312 ± 
10(syst.), using a Glauber calculation relating N part to 
the fractional cross section observed in bins of the paddle 
signal @||. 

As the geometrical layout of the PHOBOS detector 
leads to an asymmetry in the acceptance and detection 
efficiency for positively and negatively charged particles 
for a given magnet polarity, data was taken using both 
magnet polarities (BPLUS and BMINUS). The repro- 
ducibility of the absolute field strength was found to be 
better than 1%, based on Hall probe measurements for 
each polarity and the comparison of mass distributions 
for identified particles for the two polarities. 

The stability of the trigger event selection was checked 
using the number of reconstructed straight line particle 
tracks outside the magnetic field. The straight track 
multiplicities for the BPLUS and BMINUS data sets 
agreed to within 0.2%. 
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FIG. 1. Distribution of average energy loss as a function of 
reconstructed particle momentum. Three clear bands can be 
seen, corresponding to pions, kaons and protons. The solid 
lines indicate the cut regions for counting identified particles. 

To optimize the precision of the vertex and track 
finding only events with a reconstructed primary vertex 
position between -16 cm < z vtx < 10 cm along the beam 
axis were selected. We also restricted the analysis to 
tracks in the central region of the spectrometer planes 
along the 45° axis of the experiment. A precise knowl- 
edge of the primary vertex for each event is essential for 
this analysis, as the distance of closest approach of each 
reconstructed track with respect to the primary vertex 
(dca v t x ) is the only tool for rejecting background particles 
from decays and secondary interactions. By requiring 
a consistent vertex position from the SPEC and VTX 
subdetectors, in combination with the known position of 
the beam orbit, a vertex resolution better than 0.3 mm 
(RMS) in y and z directions and better than 0.5 mm in 
x direction was achieved. 

Particle tracks in the spectrometer were found as 



follows: tracks in the first 6 layers outside the magnetic 
field were reconstructed by a road-following algorithm. 




FIG. 2. Acceptance of the spectrometer as a function 
of transverse momentum and rapidity for pions, kaons and 
protons. The acceptance is averaged over the accepted vertex 
range and azimuthal angle <j). 

Inside the field, two-hit combinations on consecutive 
layers were mapped into (1/p, O) space, where O is 
the polar angle at the primary interaction vertex and 
p the total momentum. A clustering algorithm then 
selected combinations of matching hits in (1/p, O) space, 
yielding track pieces inside the field. For track pieces 
inside and outside the field, the specific ionization dE/dx 
was calculated using the truncated mean of the angle- 
corrected hit energies, discarding the two highest energy 
values. Straight and curved track pieces were then 
matched based on O, a fit to the combined track in the 
j/z-plane and requiring consistency in the independently 
determined dE/dx values. 

For the combined tracks the dE/dx was again calcu- 
lated, discarding the 4 highest energy values. Particle 
identification was performed using p and dE/dx, which 
depends only on particle velocity. The identification cuts 
are shown in Fig. [l] with three bands corresponding to 
pions, kaons and protons. The corresponding acceptance 
regions for identified particles in transverse momentum 
Pt and rapidity are shown in Fig. |^. Table | gives 
the resulting event and identified particle statistics for 
all combinations of magnet polarity, particle charge 
and particle mass, as well as the uncorrected average 
transverse momentum (pr)- Using two magnetic field 
settings allows the determination of two statistically 
independent values for each of the three particle ratios. 
As can be calculated from the numbers in Table |, the 
two values agree within statistical uncertainty for each of 
the ratios. 

Preliminary PHOBOS data show a total multiplicity 
of produced charged particles N c h ~ 3500 [[l7j for the 
event selection used here. This is large compared to the 
average initial charge asymmetry given by the number 
of participating protons of {N^ art ) w 125. The charge 
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asymmetry 1 — (N~) /(N + ) in the final state will therefore 
be small, so that a random contamination of the iden- 
tified particle samples would move the observed ratios 
towards unity. The contamination of the proton and 
antiproton samples was checked by testing the stability 
of the antiproton/proton ratio against variation of the 
particle-ID cuts shown in Fig. [j]. Within statistical error 
the ratios were stable against further changes of the cuts. 

For comparison with data from other experiments 
and theoretical calculations several corrections to the 
observed particle ratios have to be applied. These cor- 
rections account for particles produced in secondary in- 
teractions, loss of particles due to absorption in detector 
material, and feeddown particles from weak decays. Both 
feeddown and secondary particles tend to be produced 
with lower pt than the corresponding primary particles. 
Simulations show that the difference in contamination of 
7r + and 7T~ by positrons and electrons, which cannot be 
rejected within the resolution of the dE/dx measurement, 
is negligible. 

The acceptance for secondary and feeddown particles 
is limited to those produced within 10 cm radial distance 
from the primary collision vertex, as accepted tracks were 
required to have at least one hit in the first two layers 
of the spectrometer. The background contamination is 
further reduced by requiring the particle tracks to have 
dca v tx < 3.5 mm. 

As HIJING reproduces the total charged particle mul- 
tiplicity near mid-rapidity to within 10%, the simulation 
should give an accurate estimate of the background from 
secondary interactions. The correction factors are shown 
in Fig. |^ as a function of transverse momentum. The 
net correction is 0.01 for (p)/(p). For pions and kaons 
the contribution from secondary particles is much less 
than 1%. The difference in absorption for protons and 
antiprotons in the detector material was studied based on 
GEANT simulations. The fraction of absorbed protons 
and antiprotons as a function of pr is also shown in 
Fig. ||. From this correction, which averages to 1.7% for 
the protons and 6.3% for the antiprotons, a correction 
to (p)/(p) of 0.04 is obtained. The contribution from 
feeddown particles can not be modelled as precisely, as 
the absolute yield of strange hadrons has not yet been 
measured at RHIC. 



TABLE I. Number of accepted events, uncorrected number 
of identified particles and uncorrected (pr) for each magnetic 
field polarity. 
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26509 events 




41850 events 
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254 ± 1 
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FIG. 3. Correction factors for particle multiplicities as 
a function of transverse momentum px for positive and 
negative particles. Correction factors of less than 1% and 
the pr-independent feeddown correction for {p)/{p) are not 
shown. 

Simulations show that the feeddown correction for 

(K-)/(K+) and (OA""") ( see Fi S- |) is less than !%, 
whereas feeddown to protons and antiprotons from A 
and A is not negligible. The feeddown correction to 
the (p)/(p) ratio can be evaluated as a function of the 
(A)/(A) and (A) / (p) ratios. Transport model studies, as 
well as quark counting arguments and data from lower 
energies show that to a good approximation the following 
relationship holds |ll| : 

<A) = (A^+) <p> 
(A) (K-) (p) 

Our data show (K + ) / (K~) = 1.1, suggesting only a 10% 
difference in (A) / (A) compared to (p) / (p) . Using the 
dca v t x < 3.5 mm cut, less than half of the protons from 
weak decays are accepted in the analysis. This further 
reduces the feeddown correction to (p) / (p) . Using MC 
simulations, the range of the correction was estimated by 
varying the (A)/(p) ratio from 0.2, the value predicted by 
HIJING, to 0.4 and varying {K + )/{K~) from 1 to 1.2. 
The resulting correction to (p) / (p) ranges from to -0.03, 
with a most probable value of -0.01. 

After corrections we find the following ratios within 
our acceptance: 

(n-) /(tt+) = 1.00 ± 0.01(stat.) ± 0.02(syst.) 
(K-) /(K+)= 0.91 ± 0.07(stat.) ± 0.06(syst.) 
(p)/(p) = 0.60 ± 0.04(stat.) ± 0.06(syst.) 

For all three particle species the average transverse 
momentum of antiparticles and particles agrees within 
statistical error (see Table , indicating no strong depen- 
dence of the ratios on transverse momentum. The results 
are in good agreement with data recently presented by 
the BRAHMS, PHENIX and STAR collaborations [uj. 
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FIG. 4. (K~) / (K + ) and {p)/{p} ratios as a function of i/s 
for nucleus-nucleus collisions, in comparison with predictions 
from the HIJING and RQMD models. Only statistical errors 
are shown. 



using a statistical model calculation [^3| shown in Fig. 0. 
For a realistic range of freeze-out temperatures of 160 to 
170 MeV, both (K + )/{K~) and (p)/(p) are consistent 
with fiB = 45 ± 5 MeV. This is much lower than the 
value of fi B = 240 - 270 MeV [|J|] obtained in statistical 
model fits to Pb+Pb data at ^/s~^ = 17.2 GeV, showing 
a closer but not yet complete approach to a baryon-free 
regime at RHIC. 
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In Fig. |4| our results are compared to lower energy data 
3 |H and calculations using the HIJING and RQMD 
| microscopic transport models. The values for the 
I<~)/{K + ) and (p)/(p) ratios are significantly higher 
than at lower energies. Using the same parameters that 
successfully described the charged particle multiplicity 
density in Au+Au collisions at RHIC 0], HIJING 
overestimates (p)/(p) by 0.18. Clearly, the particle 
ratio data provide additional constraints for microscopic 
models. A better description may require a modification 
in the baryon number stopping mechanism, baryon-pair 
production or antibaryon annihilation. An indication 
of possible mechanisms is given by the comparison with 
RQMD, which includes rescattering of produced hadrons 
and predicts a value for (p)/{p) that is 0.18 below our 
data. 
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FIG. 5. Statistical model calculation (dotted lines) of 
(K~) / (K + ) and (p}/(p) as a function of fiB by Redlich et 
al. The horizontal bands show the ratios observed in the data 
(statistical errors only). The vertical shaded area indicates 
the allowed region in \ib- 

Finally, we estimate the baryo-chemical potential [Ib 



[1] See e. g. J. P. Blaizot, Nucl. Phys. A661 (1999) 3. 

[2] E. Schnedermann, J. Sollfrank, U. Heinz, Phys. Rev. C48 

(1993) 2462. 
[3] F. Becattini, Z. Phys. C69 (1996) 485. 
[4] P. Braun-Munzinger, I. Heppe, J. Stachel, Phys. Lett. 

B465 (1999) 15. 
[5] P. Koch, B. Miiller, J. Rafelski, Phys. Rep. C142 (1986) 

167. 

[6] S. A. Bass et al., Prog. Part. Nucl. Phys. 41 (1998) 225. 
[7] H. Sorge, Phys. Rev. C52 (1995) 3291. 
[8] W. Busza, A. S. Goldhaber, Phys. Lett. B139 (1984) 
235. 

[9] D. Kharzeev, Phys. Lett. B378 (1996) 238. 
[10] S. E. Vance, M. Gyulassy, X. N. Wang, Phys. Lett. B443 
(1998) 45. 

[11] B. Back et al., Nucl. Phys. A661 (1999) 690. 
[12] H. Pernegger et al., Nucl. Instrum. Methods. A419 
(1998) 549. 

[13] R. Pak, PHOBOS Collaboration, proceedings of Quark 

Matter '01 Conference (in press). 
[14] B. Back et al., Phys. Rev. Lett. 85 (2000) 3100. 
[15] J. Katzy et al., PHOBOS Collaboration, proceedings of 

Quark Matter '01 Conference (in press). 
[16] M. Gyulassy and X. N. Wang, Comp. Phys. Comm. 

83 (1994) 307. We used HIJING V1.35 with standard 

parameter settings. 
[17] A. Wuosmaa et al., PHOBOS Collaboration, proceedings 

of Quark Matter '01 Conference (in press). 
[18] J. Zimanyi, Nucl. Phys. A661 (1999) 224. 
[19] See proceedings of Quark Matter '01 Conference (in 

press) . 

[20] L. Ahle et al., Phys. Lett. B490 (2000) 53, 

Phys. Rev. C60, (1999) 064901, 

Phys. Rev. Lett. 81, (1998) 2650 . 
[21] J. Bachler et al., Nucl. Phys. A661 (1999) 45. 
[22] I. G. Bearden et al., Phys. Lett. B388 (1996) 431. 
[23] K. Redlich, proceedings of Quark Matter '01 Conference 

(in press). 



4 



